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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
in  Partial  Fulfillment  of  the  Requirements  for  the  Degree  of 

Doctor  of  Philosophy 

DYNAMIC  BEHAVIOR  OF  HEATED 
CYLINDRICAL  AND  CONICAL  SHELL§ 

By 

Chung-li  Sun 
August,  1967 

Chairman:  Dr.  S.  Y.  Lu 

Major  Department:  Engineering  Science  and  Mechanics 

The  present  analysis  is  to  study  the  elastic,  dynamic  stability 
and  response  of  conical  and  cylindrical  shells  subjected  to  time- 
dependent  heating.  Two  phases  of  the  problem  were  treated.  First, 
a linear  shallow  shell  theory  was  used.  The  shells  become  dynamically 
unstable  when  the  normal  displacement  component  due  to  thermal  stresses 
increases  unboundedly  with  time.  Second,  a nonlinear  shallow  shell 
theory  has  been  adopted.  The  post-buckling  amplitude-frequency  relation- 
ship shows  that  the  solution  originating  at  the  upper  boundary  of  the 
instability  region  is  always  unstable,  while  the  solution  originating 
at  the  lower  boundary  of  the  instability  region  is  always  stable. 

Two  kinds  of  temperature  change  with  respect  to  time  were  con- 
sidered. In  one  case,  the  temperature  change  is  sinusoidal,  and  in 
the  other  case,  the  temperature  function  is  a ramp-type.  The  shells 
are  simply  supported  and  completely  restricted  at  both  ends.  The 
Hamilton  principle  has  been  used  for  both  the  linear  and  nonlinear 
cases.  The  final  equation  in  both  cases  is  of  the  Mathleu  type.  The 
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linear  Mathleu  differential  equation  was  used  to  study  the  charac- 


teristics of  the  dynamic  stability  of  the  shells,  and  the  nonlinear 
Mathleu  differential  equation  was  used  to  establish,  the  amplitude- 
frequency  diagram,  so  that  the  stable  and  unstable  solutions  were 
thus  determined. 
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CHAPTER  1 


INTRODUCTION 

In  recent  years,  owing  to  the  rapid  development  of  structural 
engineering,  a new  field  entitled  dynamic  stability  of  elastic  plates 
and  shells  has  been  intensively  cultivated.  The  theory  used  in  this 
field  is  generated  from  the  study  of  stability  of  elastic  systems 
under  loads  varying  with  time.  In  the  analysis  of  dynamic  stability 
of  deformable  bodies,  the  influences  of  inertia  must  be  taken  into 
consideration,  and  the  system  becomes  dynamically  unstable  when  the 
amplitudes  of  the  vibrations  increase  unboundedly  with  time.  Thus, 
the  theory  of  dynamic  stability  of  elastic  systems  can  be  defined  as 
the  study  of  vibrations  induced  by  pulsating  or  impact  loadings  which 
are  parametric  with  respect  to  certain  displacement  components. 

This  point  of  view  can  be  illustrated  by  considering  a rod  com- 
pressed by  a time-dependent  force.  If  a straight  rod  is  subjected  to 
a periodic  longitudinal  load  and  if  the  amplitude  of  the  load  is  less 
than  that  of  the  static  buckling  value,  then,  in  general,  the  rod 
experiences  only  longitudinal  vibrations.  However,  it  has  been  known 
that  for  certain  relationships  between  the  force  frequency  and  the 
natural  frequency  of  transverse  vibration  a straight  rod  becomes  un- 
stable, and  transverse  vibrations  occur.  The  amplitude  of  this 
vibration  increases  rapidly  to  large  values . Such  a phenomenon  is 
called  parametric  resonance.  The  relationship  of  the  frequencies  at 
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such  a resonance  differs  from  the  frequency  relationship  for  the 
ordinary  forced  vibration  resonance.  For  sufficiently  small  values 
of  the  longitudinal  force,  this  relationship  is  that  the  force  frequency  is 
equal  to  twice  the  value  of  natural  frequency  of  the  rod. 

Similarly,  periodic  forces  acting  in  the  middle  plane  of  a thin 
plate  or  shell  can  excite  intensive  transverse  vibrations  under  certain 
conditions. 

Another  example  is  when  a periodic  loading  is  applied  symmetrically 
with  respect  to  an  arch.  It  causes  not  only  symmetrical  vibrations,  but 
sometimes  can  excite  asymmetrical  vibrations  of  very  large  amplitude. 

In  general,  whenever  there  is  a special  kind  of  static  loading 
which  causes  a loss  of  static  stability,  the  vibrational  loading  of  the 
same  type  will  cause  a loss  of  dynamic  stability.  Such  a loading  is 
characterized  by  the  fact  that  it  is  contained  as  a parameter  in  the 
equation  of  motion.  Therefore,  it  is  called  a parametric  loading.  This 
term  is  more  appropriate  because  it  indicates  the  relation  to  the  phe- 
nomenon of  parametric  resonance. 

The  general  concept  stated  above  has  been  discussed  by  Ziegler  [1] 
and  Bolotin  [2].  In  the  latter  reference,  the  general  theory  of  the 
dynamic  stability  of  elastic  systems  is  formulated,  and  examples  for 
the  cases  of  rods,  beams,  frames  and  plates  under  the  periodic  loading 
have  been  given.  Linear  problems  of  the  circular  cylindrical  shell 
under  a uniformly  distributed  radial  periodical  loading  as  well  as  a 
compressive  periodical  loading,  and  spherical  shells  subjected  to  a 
radial  periodical  loading  uniformly  distributed  over  the  surface  have 
been  discussed  also  in  reference  [2], 

A correlated  problem  to  the  problem  of  dynamic  stability  is 


3 


the  dynamic  buckling.  An  elastic  structure  under  any  time-dependent 
loading  may  cause  a finite  Jump  of  the  deflection  at  certain  value 
of  the  loading.  This  value  is  defined  as  the  dynamic  buckling  load. 

For  structures  with  large  initial  imperfections,  such  a jump  will  not 
happen.  However,  there  would  be  a point  of  inflection  on  the 
deflection-loading  curve,  and  the  value  of  the  loading  at  this  point 
is  considered  to  be  the  dynamic  buckling  load.  A general  criterion 
for  dynamic  buckling  of  elastic  structures  which  are  sensitive  to 
initial  imperfections,  and  some  estimates  for  the  dynamic  buckling 
strengths  of  such  structures  subjected  to  various  loading  histories 
can  be  found  in  the  paper  of  Budiansky  [3]. 

The  study  of  the  problem  of  dynamic  stability  began  in  the 
early  days  of  this  century.  In  1924,  Beliaev  [4]  examined  a problem 
of  a hinged  straight  rod  under  axial  loading.  Later  in  1935,  Krylov 
and  Bogoliubov  [5]  studied  the  same  problem  with  arbitrary  support 
conditions  by  applying  the  Galerkin  method.  In  this  part  of  the 
world,  Lubkin  and  Stoker  [6]  first  investigated  this  kind  of  problem 
in  1943,  by  analyzing  the  stability  of  columns  and  strings  under 
periodically  varying  forces.  Hoff  [7],  in  1954,  examined  the  stability 
of  bars  with  initial  imperfection  under  dynamic  loading  by  a com- 
pressing force. 

In  recent  years,  several  authors  studied  the  linear  problems 
of  dynamic  stability  of  shells  in  further  detail.  Wood  and  Koval  [8], 
in  1963,  explored  some  experimental  and  analytical  results  of  the 
collapse  of  thin-walled  cylindrical  shells  with  the  initial  imper- 
fection under  dead  weight  axial  compression  and  axisymmetric  transient 
and  oscillatory  hydrostatic  pressures.  An  eighth  order  partial 
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differential  equation  of  the  Donnell  type  [9]  has  been  solved  by 
assuming  a double  Fourier  series  for  the  solution  to  the  normal 
displacement.  The  resulting  second  order  differential  equation  of 
the  coefficients  of  the  displacement  function  with  respect  to  time 
was  discussed  on  its  stability  under  the  step-pressure  or  sinusoidal 
pressure.  In  the  same  year,  Yao  [ 10 ] studied  the  dynamic  stability 
bounds  for  complete  cylindrical  shells  under  various  combinations 
of  static  and  periodic  radial  as  well  as  axial  loads.  The  solutions 
were  determined  numerically  by  means  of  properties  of  the  Mathieu 
functions.  Goodier  and  Mclvor  [ll],  in  1964,  considered  the  surface 
of  a thin-walled  cylinder  subjected  to  impulsive  pressure.  They 
assumed  the  form  of  the  normal  displacement  and  substituted  it  into 
the  formulated  Lagrange  equation  of  motion.  The  stable  and  unstable 
modes  of  the  normal  displacement  were  discussed  through  the  obtained 
Mathieu  equation.  In  1966,  Bieniek,  Fan  and  Lackman  [12]  introduced 
a perturbation  method  to  solve  the  same  problem  under  ramp-type 
pressure  and  sinusoidal  pressure.  The  result  of  the  ramp-type 
pressure  agrees  with  that  in  reference  [ll]. 

A nonlinear  problem  of  dynamic  stability  of  spherical  shell  was 
investigated  in  1958  by  Bolotin  [13].  In  this  reference,  the  un- 
Bymmetric  character  of  the  nonlinearity  on  a curved  panel  reveals  the 
essential  difference  in  the  behavior  of  a plate  and  a curved  panel 
for  parametric  excitation.  In  the  same  year,  Vol'mir  [14]  by  using 
the  Galerkin  method  and  the  method  of  numerical  integration  solved 
the  problem  of  dynamic  stability  of  a cylindrical  panel  with  initial 
imperfection  subjected  to  axial  compression.  In  1959,  Agamirov  and 
Vol'mir  [15]  used  the  same  approximate  technique  as  used  in  reference 
[14]  to  treat  the  same  problem  of  a complete  thin  cylinder  under  a 
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linearly  increased  time-dependent  hydrostatic  pressure.  In  both 
papers,  the  solution  to  the  transverse  displacement  is  assumed  to 
the  first  degree  of  approximation,  and  a second  order  nonlinear 
differential  equation  of  the  time-dependent  amplitude  of  the 
normal  displacement  was  thus  obtained.  After  this  differential 
equation  was  numerically  integrated,  the  critical  conditions  were 
depicted  by  the  characteristics  of  the  response  curve.  In  1961, 
Mishenkov  [16]  considered  a shallow  cylindrical  panel  under  axial 
periodical  compression.  The  Galerkin  method  wsb  applied  again  by 
representing  the  normal  displacement  by  the  first  term  of  a double 
Fourier  series  with  time-dependent  coefficient.  A nonlinear  differ- 
ential equation  of  the  second  order  with  periodic  coefficient  was 
finally  obtained,  and  periodic  solutions  to  this  equation  near 
principal  resonance  were  solved.  The  stability  was  thereupon  in- 
vestigated. In  1962,  Cummings  [17]  examined  both  the  cylindrical 
panel  and  the  complete  thin  cylindrical  shell  under  different  forms 
of  radial  pressure.  Although  the  investigation  was  interested  mainly 
in  the  problem  of  vibration  and  response,  the  characteristics  of  loss 
of  stability  in  a linear  sense  were  also  mentioned.!  In  1964,  Ho  and 
Nash  [18]  discussed  the  dynamic  instability  of  thin  shallow  spherical 
shells  under  uniformly  distributed  impulsive  loadings  including  an 
impulsive  function,  a step-function  and  a linear  function  of  time. 
They  used  the  Galerkin  method  as  well  as  a dynamic  stability  theory 
to  establish  their  solutions.  The  same  techniques  were  then  extended 
to  the  solution  of  complete  spherical  shells.  In  1965,  Wilder  and 
Nash  [19]  investigated  the  response  of  a truncated  conical  shell 
subjected  to  the  dynamic  axial  force.  After  first  applying  the 


6 


variational  principle  to  find  out  the  nonlinear  differential  equation 
of  the  second  order  with  respect  to  time,  the  Runge-Kutta  numerical 
integration  method  was  adopted  to  trace  out  the  nonlinear  response. 

All  papers  reviewed  above  followed  the  important  simplification 
suggested  by  Reissner  [20]  in  1955,  that  the  tangential  inertia  terms 
for  thin-walled  surface  structures  may  be  omitted  with  negligible 
errors . 

Although  the  study  of  dynamic  stability  of  cylindrical  shells, 
by  using  both  linear  and  nonlinear  shallow  shell  theories,  has  been 
developed  in  such  a broad  scope,  the  parallel  study  of  a conical  shell 
still  remains  in  a comparatively  primitive  stage.  Moreover,  similar 
problems  dealing  with  a temperature  loading  are  even  less  discussed. 

Only  a few  problems  of  the  linear  dynamic  response,  and  certain  problems 
of  the  static  stability  of  beams,  plates  and  shells  under  temperature 
loadings  are  listed.  Boley  [21]  first  treated  the  problem  of  thermally 
induced  vibrations  of  beams  in  1956.  In  this  paper,  he  introduced  a 
thermal  moment  in  the  beam  equation  by  applying  a sudden  heat  input  to 
the  upper  surface  of  the  beam  uniformly  distributed  along  its  span.  The 
ratio  of  the  dynamic  response,  obtained  by  solving  this  equation,  and 
the  corresponding  quasistatic  deflection  was  plotted  against  a governing 
nondimens ional  time-ratio  parameter.  From  this  curve,  Boley  concluded 
that  the  maximum  deflection  under  the  dynamic  effect  is  twice  larger  than 
that  calculated  through  a quasistatic  process.  In  1957,  Boley  and 
Barber  [22]  extended  the  same  problem  to  a plate,  and  drew  out  an  analogous 
conclusion.  Later  in  1966,  Kraus  [23]  investigated  the  same  problem  for 
spherical  shells  by  use  of  a model  expansion  procedure.  In  1967,  Lu 
and  Sun  [24]  examined  this  problem  for  truncated  conical  shells  through 
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a variational  method.  The  significance  of  dynamic  effects  on  deflection 
are  evident  in  these  references  [21-24]. 

Hoff  [25],  in  1957,  solved  the  thermal  buckling  problem  of  hinged 
thin  cylindrical  shells  subjected  to  uniform  surface  temperature  loading. 
He  used  the  Donnell  equation  and  concluded  that  the  elastic  thermal 
buckling  is  unlikely  to  occur  except  when  the  shell  is  very  thin.  In 
1962,  Sunakawa  [26]  studied  the  same  problem  for  a clamped  thin  cylinder 
and  found  that  the  shell  will  buckle  at  a lower  temperature  than  that 
for  simply  supported  ends.  In  1959,  Huang  and  Van  Der  Maas  [27]  worked 
out  the  buckling  problem  of  heated  and  simply  supported  plates  with  the 
Von  Kirmfinn  large  deflection  theory  and  Rayleigh-Ritz ' s method.  Solution 
forms  of  the  three  components  of  displacement  were  chosen  to  be  the  first 
term  of  a double  trigonometric  series.  For  the  two  inplane  displacement 
components,  linear  terms  were  added  to  justify  the  boundary  conditions. 

In  1960,  Boley  [28]  wrote  a book  on  thermal  stress  analysis,  in  which 
both  basic  theories  and  applications  to  rods,  beams,  plates  and  shells 
are  presented.  In  1964,  Hoff,  Chao  and  Madsen  [29]  used  the  Donnell 
equation  to  analyze  the  buckling  problem  of  a thin-walled  cylindrical 
shell  heated  along  an  axial  strip.  In  this  report,  solutions  for  heated 
and  unheated  regions  were  obtained  separately,  then  coupled  together 
through  continuity  of  boundary  conditions.  Both  symmetric  and  anti- 
symmetric bucklings  were  discussed.  The  conclusion  is  that  the  critical 
value  of  the  compressive  axial  stress  is  found  to  be  equal  to  the 
critical  stress  of  the  same  circular  cylindrical  shell  when  subjected 
to  uniform  axial  compression  if  the  band  is  not  extremely  narrow.  Other- 
wise, the  critical  stress  will  be  higher. 

The  proposition  of  this  dissertation  consists  of  two  parts.  In 
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the  first  part,  a linear  theory  of  shallow  shells  is  adopted.  Some 
qualitative  characteristics  of  the  dynamic  stability  of  thin  cylindri- 
cal and  conical  shells  related  to  the  induced  time-dependent  thermal 
loading,  the  geometry  of  the  shells  and  the  mode  number  of  the  normal 
displacement  are  investigated.  In  the  second  part,  a nonlinear  theory 
of  shallow  shells  is  followed.  The  post-buckling  thermal  response- 
frequency  relationship  was  examined  so  that  the  stable  and  unstable 
region  could  be  defined  through  the  "jump"  phenomenon. 

The  shells  under  consideration  are  assumed  thin  and  of  constant 
thickness.  Two  kinds  of  temperature  change  with  respect  to  time  were 
considered.  In  one  case,  the  temperature  gradient  is  sinusoidal,  while 
in  the  other  case,  the  temperature  function  is  a ramp-type.  The  shells 
are  simply  supported  and  completely  restricted  at  both  ends.  The 
Hamilton  principle  was  applied  to  obtain  the  differential  equations  for 
the  generalized  coordinates.  The  problem  was  treated  in  two  steps,  of 
which  the  first  is  prior  to  buckling  and  the  second  is  during  buckling 
for  the  linear  case  and  after  buckling  for  the  nonlinear  case.  The 
prebuckling  relations  are  used  to  find  the  thermal  stresses.  The  final 
differential  equation  for  the  linear  case  assumes  the  Mathleu  type,  and 
for  the  nonlinear  case  assumes  the  type  of  a second  order  nonlinear 
differential  equation  with  variable  coefficient.  This  nonlinear  differ- 
ential equation  is  of  the  same  form  as  that  of  a nonlinear  dynamic 
stability  problem  of  a plate  under  axial  periodic  loadings  [2],  The 
quadratic  nonlinearity  which  shows  up  in  the  final  differential  equation 
of  a curved  panel  is  missing  in  the  present  case  because  of  the  complete- 
ness of  the  geometry  of  the  shell  along  the  circumferential  direction. 

In  the  linear  case,  the  region  of  stability  at  various  normal 
displacement  modes  has  been  studied.  The  behavior  of  stability  depends 
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upon  the  temperature  function  as  well  as  the  geometry  of  shells.  In 
the  nonlinear  case,  the  force  is  always  a type  of  hard  spring  force  [30] 
for  the  present  boundary  conditions.  The  solutions  originating  at  the 
lower  boundary  are  therefore  always  stable  while  those  originating  at 
the  upper  boundary  are  always  unstable. 


CHAPTER  2 

BASIC  RELATIONS  AND  METHOD  OF  APPROACH 
2.1  Coordinate  System 

The  coordinate  system  chosen  for  a conical  shell  on  its  middle 
surface  is  (*,<£)  as  shown  in  Figure  1.  This  coordinate  system 
assumes  the  geometric  orthogonality  of  the  shell.  The  transformations 
between  this  set  of  coordinates  and  the  Cartesian  coordinates  (X.Y,  Z ) 
according  to  Figure  1 are  given  as  follows 


x = * Sinecoi^  ( * Sin©Sin<p,  Z=*co5© 


(2.1) 


Substituting  the  above  expressions  into  the  relation 


and  calculating  the  values  of 


. z _ cyx  t Yx 
A ” a*  bx  f 


p 2 _ 3*  2>x 
D d*> '3^ 


one  can  immediately  obtain  the  scale  factors 


A = I, 


B = ?<  Si'nQ 


(2.2) 


For  the  limiting  case  as  00  and  ©40 , i.e.,  for  the  case  of  circular 
cylindrical  shell 


A = I , B=Rl  = Ri  = R 


(2.3) 
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Also  from  Figure  1,  one  can  easily  observe  that  the  principal 
curvature  along  X and  <p  directions  are 


= ( 9C  fan  0 ) ' (2.4) 

For  a cylindrical  shell 

^X=  0,  = R (2.5) 


2.2  Elastic  Stress-Strain-Displacement  and  Moment-Curvature- 
Displacement  Relations 

In  the  present  investigation,  the  following  basic  assumptions 
are  made 

(1)  Elements  normal  to  the  unstrained  middle  surface  of 
the  shell  remain  normal  to  the  strained  middle  surface. 

(2)  The  material  obeys  Hooke's  law. 

(3)  In  the  nonlinear  analysis,  only  normal  displacement 
is  assumed  large. 

Based  upon  the  above  assumptions,  the  strain-displacement 
relations  to  the  middle  surface  of  the  shell  under  consideration 
are  [31] 


du  , I / 3w\2 

= ^ + ziTt) 


(2.6a) 


(2.6b) 


f.  r j>B.  , I 

2 l d*  ‘ & d<p  B d*  B df> 


(2.6c) 
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The  stress-strain-displacement  relations  are  [31] 


J=_  ( £ -l - d a-  ) — 


dv  , u.  £B 
B dp  B d* 


l / i 3w\2 


E(XT 
1 -U 


E<XT 

I 


E f i 3 v , u__aB_  , , j_  r j_  q\n  \ 

I B dp  + B <3;x  2 ^ B <3<£' 


I / I \2-  j 


+ 


E°(T 
I -y 


| -+  y 

_ E / dv  _l_  v_  <^E>  + j_  cM  3w  n 

“ 2(i+f) ' dx  & d<p  0 £2  B 

The  curvature-displacement  relations  are 


A/. 

d*z 


C !»-'  = 


J_  a w 

b2  a^>2 


j_  <)B  3 w 
B d-x  dx 


l d w <)  B _ _l_  buVJ 

5Z  dp  a*  b a*a<£ 


(2.7a) 


(2.7b) 


(2.7c) 


(2.8a) 


(2.8b) 


Crx^)"'=  C r^x)"' 


(2.8c) 
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The  moment-curvature-displacement  relations  are 


m*  = jd^  pc  t>r'j 


dZw  I b \aJ 

T*  "TT* 


d *2  ^ &2  d<pz  B dx  d* 


)) 


(2.9a) 


M<ji=  d(  ( W w r^r'3 


_ «•_[_  dzw  , J_  d&  ^ w d2  w 
“ B S'*2- 


) 


(2.9b) 


M^=  =D(I-W(r^) 


-I 


= D0-jO( 


J_  b w gB 

B2  b<p 


j_  ()2  w 'I 

B 


(2.9c) 


In  the  above  expressions,  u,  v,  w are  the  displacement  components 

at  an  arbitrary  point  on  the  middle  surface  of  the  shell  and  are 

measured  along  %<£  and  normal  directions,  respectively  (Fig.  1). 

The  notations  £*,  6^  , 6^  and  <rx , , cTy<^  are  the  normal  and  shearing 

strains  and  stresses  to  the  middle  surface;  M^,M^  and  are  the 

bending  and  twisting  moments,  while  C rx)  ' ( fy) ' (rx^)  represent  the  changes 

of  curvature.  In  this  dissertation,  B assumes  Young's  modulus,  V is 

Poisson's  ratio,  cX  represents  the  linear  thermal  expansional  coefficient 

and  T=T(  X , <p , X-j , i ) is  a temperature  function  in  which  Xj  is  the  co- 

Eh* 

ordinate  along  the  normal  direction,  and  t is  time.  D- 
which  h denotes  the  thickness  of  the  shell. 
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2.3  Energy  Expressions 

Referring  to  the  middle  surface  of  the  shell,  the  strain  energy  u 
can  be  considered  as  the  sum  of  three  parts  [31],  i.e., 

u - 

where  17^ ; Ub , rePreaent  the  membrane,  bending  and  thermal  energy, 
respectively.  The  expressions  of  these  parts  of  energy  in  terms  of 
strains  and  changes  of  curvature  and  finally  in  terms  of  displacements 
can  be  written  as  follows  [31] 


= + + + Z (i'P)e%<pJ  Bdxdcp 

s 

~ 2 ^ +-  -g  w + 

5 

+ 2 w^)2J2+2i>t  U,x  + 75 


w,4,)2J  + j( 1 vn  + 


l . -t  2 

+ "g  W;<^J 


B d9f 


(2.10) 


Vb - £ JJ  C ^ l)  + (•  r^1)  + 2.  v ( r^')(r^')n-2 ( i-u)( rx$)2]  0 c|t< d^> 


= ^Jf  t ( + ( ~q  B>>  w>*  )2  + 

5 

+ 2plW,rt(“j  + 

-t  B W>*£_  3 6 d*d^>  (2.11) 
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= 


■ ■“  jj"  t ( ex+^4!)T  -*-  ( r'1  + r^‘)  mt]  B d-* 


i -yJJ  ' [ ^;3f  + 2 " w’*'  + b ^ ^ 3®;x  w T 


p ( p,  W,^)  J T "*■  (.  '*’ 


-t 


"g  mt}  b 


(2.12) 


where  the  integral  is  over  the  total  middle  surface  S • The  notations 


T and  M -r  are  defined  by  the  integrals 


T 


=ih  *Td 


X-. 


Mr  = /h  *j*T<for3 


The  symbol  D'  used  above  is  defined  by  15  = ~jj2  • The  subscript  ( , ) 
represents  the  partial  differentiation  with  respect  to  the  following 
coordinates. 

In  case  T is  not  a function  of  #3 


Vr  =-777//  (<?x+e<t> ) t 6d*d<f 
s 

= --p£.f/  t +yC^n)2+  eV'tB>*-,Vv  + 

5 

+ w),4,)2j«TB 

The  kinetic  energy  of  the  shell  can  be  expressed  as 

V = y //  f //  C ( “ + *3  )Z  + ( V/  + *3  ^3} 


(2.13) 


=f//i 


h(U  + V +w 


B cU  d 4> 

• 2 • 2 , , li  \z.  h3  / ^;4M2 


J + 1T^/+  &(*£■) 
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where  the  dot  ( * ) represents  , f is  the  mass  density.  Since  the 
shell  is  assumed  thin,  the  rotatory  inertia  terms  can  be  reasonably 
neglected.  Also  for  thin  shallow  shells,  as  Reissner  [20]  pointed 
out  in  his  work  that  the  tangential  inertia  forces  can  be  ignored 
with  negligible  errors,  the  inertia  effect  considered  in  this 
dissertation  is  limited  only  to  the  inertia  effect  along  the  normal 
direction.  Thus, 


V = ^ //  C w)z  B cl*  d4>  (2.15) 

S 

For  solving  the  present  shell  problems,  the  Hamilton  principle 
is  adopted.  This  principle  states  that  among  all  the  admissible 
states  of  stress  the  actual  one  is  that  which  satisfies 


(V-U)dt  = o 


(2.16') 


where  8 is  a variational  symbol,  and  U,  V are  the  strain  energy  and 
the  kinetic  energy,  respectively. 

Expression  (2.16')  can  be  arranged  to  a more  suitable  form 
for  applying  the  direct  variational  method.  Since 

8 j*1  Vdt  = 8 [tz(  If— ■f*)z5d'xd4]dt 

[ J / fh  w Sw  6 dxd(f>]  dt 

r,  S 

One  can  integrate  this  expression  by  parts  with  respect  to  t,  to  yield 


<s/t!  Vd+ = r\  ji-fh 


( i a/  <5w)  + 


^ ( Jf~  fh  ( kv  <5  vj)  6 dxd<p  J 2 

* i 
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The  variation  of  the  normal  displacement,  for  the  present  problem, 
is  always  assumed  zero  at  both  terminals  of  the  time  interval. 
Therefore,  it  gives 

hf  1 Vdt  = f [ J f~  f h ( iv  <§ia/)  bdKd<p]  di 

s 

Thus,  Equation  (2.16')  can  be  written  in  a new  form 

J,  [Jf  fhfwSvMjBdydcp  + 8 If]  dt  (2.16) 

' 5 

This  expression  is  the  actual  one  which  will  be  used  to  carry  out 
the  direct  variational  calculations  in  the  present  dissertation. 


2.4  Boundary  Conditions  and  Initial  Conditions 

In  this  dissertation,  the  boundary  conditions  are  simply 
supported  and  are  completely  restrained  at  both  ends.  These 
conditions  can  be  formulated  as 


_ ( X = X , X,  (conical  shell) 

w = U = V = 0 at  I ’ 2 

= t L/2  (cylindrical  shell) 

, o(  = ?(  , X?  (conical  shell) 
Mx=-D(rr,+  yr2-,)=o  at  ( 

' *x  — i *-/2.  (cylindrical 
shell) 

where  is  the  longitudinal  bending  moment. 

The  initial  conditions,  for  t =0  , are  assumed  to  be 


(2.17) 


W = w - O (2.18) 

2.5  Method  of  Approach 

As  it  has  been  stated  in  the  first  chapter,  the  solution  is 
carried  out  in  two  steps.  The  displacements,  stresses  and  moments 
in  the  shell  before  buckling  are  found  first.  In  the  prebuckling 
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state,  the  deformations  are  assumed  symmetric.  Thus,  v'=o  , and 
U and  w depend  on  (o^t  ) only.  The  superscript  ( ’ ) is  applied 


to  quantities  at  unbuckled  state.  The  second  step  of  the  solution 
is  to  study  the  buckling  relations  for  the  linear  case  and  the  post- 
buckling  relations  for  the  nonlinear  case.  The  superscript  ( *' ) is 
used  to  define  the  additional  quantities  due  to  buckling,  or  after 
buckling.  The  additional  displacement  components  are  denoted  by 
( «*  v'l  w").  They  are  the  functions  of  ( 'X,  4s,  "t  ).  Both  ( u’,  uv/') 
and  ( u",  V",  \A’  ) are  chosen  to  satisfy  the  boundary  conditions. 

It  has  been  known  that  the  shell  deformations  before  buckling 
agree  well  with  the  classical  theory  of  shells.  This  means  that 
the  quadratic  terms  of  the  strain  components  £*,£<*,,  6 *<f>  have  nothing 
but  a negligible  influence  to  the  solution  of  prebuckling  deformations, 
stresses  and  moments.  Therefore,  the  strain  components  can  be  separated 
into  two  parts  t 6*  ; tq,  and  6x(j,  ; of  which  all  the  strains  with 
the  superscript  ( ’ ) are  linear  strains,  and  those  with  the  superscript 
(’  ) are  the  additional  nonlinear  strains.  They  can  be  expressed  as 


(2.19a) 


(2.19b) 


(2.19c) 


(2.20a) 


(2.20b) 


(2.20c) 
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For  the  prebuckling  solution  of  deformations,  stresses  and  moments, 
ex'  , €-4,  and  6*4,  instead  of  €rx  , €<p,  are  substituted  in 
Equations  (2.7),  (2.10)  and  (2.13).  Then,  according  to  the  conditions 
(2.17),  the  displacement  functions  u and  w'  are  assumed.  The  ampli- 
tudes of  u’  and  w’  are  functions  of  time  and  known  as  generalized 
coordinates.  Let  u'  and  w’  replace  Uandw,  respectively,  in  the 
modified  Equations  (2.10),  (2.13)  and  (2.11),  (2.15).  After  inte- 
gration over  the  total  middle  surface,  the  Hamilton  principle  is 
applied  to  establish  the  equations  needed  for  solving  the  generalized 
coordinates  of  U.;  and  vj’  , The  thermal  stresses  ( (T/,  <?$  , (X x<p  ) and 
the  thermal  moments  ( MX  / ) are  thus  determined  from  the 

modified  Equations  (2.7)  and  (2.9),  with  V = o and  u,w  replaced  by 

a’,  w’  • 

During  buckling,  the  slope  changes  of  the  shell  deformation 
increase  rapidly,  and  additional  stresses  ( (rj  ; (T^  , ) and 

//  /'  It 

additional  moments  ( M*  , Mq>/  ) are  induced.  At  this  stage,  the 
classical  theory  of  shells  can  no  longer  handle  this  problem.  The 
additional  strain  energy  developed  from  the  influences  of  the  non- 
linear strain  effect,  and  the  additional  stresses  and  moments  must  be 
introduced  to  justify  the  result  of  the  present  problem.  Nevertheless, 
the  direct  coupling  terms  of  the  nonlinear  strains  and  the  additional 
stresses  still  remain  to  be  small  quantities  of  higher  order,  and  can 
be  neglected  at  this  stage  without  serious  error.  Thus,  one  can 
formulate  the  membrane  strain  energy  during  buckling  by  the  following 
steps.  First,  the  additional  stresses  and  nonlinear  strains  are  added 
to  the  prebuckling  stresses  and  linear  strains  to  formulate  the  membrane 
strain  energy  as 


(zz-z)fapgfy**J}(«''^+ pflciz  ) + /*?)]  JJt  4- 

J (U 

+ +p*pg  + 

+ ^p^pgi^^z  =«n 


*^XI®UTJ  puB 


. s 

4weCi!(Vi)(«->>z*^!'’<'t+z(*i)+/*3))ff2§4. 

4-  4>p lop  g (4‘?  *>JD  2 4-  *3  *0  4-  ; j-/?  H 4. 

s 

+ 9 C *V9  y-0  £ 4-  ^9  *0  -I-  *9  *J)  J jj  £q  =Wilg 
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As  far  as  the  bending  strain  energy  is  concerned,  one  can  formulate 
it  in  the  following  straightforward  form 

8vb  = ff  (7  M*  -t  M"  ) &(r~')  + ( M+)  $(r$)  + 

s 

-f  2 ( Mx<f>  + Mx<p)d(rxQ  )J  6 ciycitp 

in  which  ( ) are  the  prebuckling  moments.  If  the 

variational  symbol  <5  is  removed  as  it  has  been  done  in  deriving  the 
membrane  strain  energy  expression  during  buckling,  the  bending  strain 
energy  will  assume  the  following  form 

~^b  ~ SJ  l r*  1 + Md>  r4>  + 2 Mx4>  rx<p  ] + 

5 

3 (2.23) 

The  displacement  components  u"  v",  u\ j"  are  assumed,  and  their 
amplitudes  are  again  functions  of  time.  Substituting  the  assumed 
u'  \j'i  w"  in  Equations  (2.8),  (2.19),  (2.20)  and  (2.15),  then  inte- 
grating  Equations  (2.15),  (2.22),  (2.23)  over  the  total  middle  surface 
of  the  shell,  the  strain  energy  during  buckling,  , and  the 

kinetic  energy  during  buckling,  V , will  be  obtained.  Applying  the 
Hamilton  principle  to  the  new  U and  new  ~\J  , a second  order  linear 
differential  equation  for  the  generalized  coordinate  of  iv  with 
variable  coefficient  is  obtained.  This  differential  equation  will  be 
of  the  Mathieu  type  in  this  dissertation.  So  the  dynamic  stability 
of  the  shells  calculated  through  a linear  theory  of  the  shallow  thin 
shell  will  be  discussed  according  to  its  properties  shown  in  the 
Mathieu  diagram. 
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After  the  buckling,  the  slope  changes  of  the  shell  deformation 
become  larger  and  larger;  also,  the  additional  stresses  and  moments 
tend  to  be  larger  and  larger.  In  this  case,  large  deformation  theory 
(nonlinear  theory)  of  shallow  thin  shells  must  be  used  to  establish 
the  strain  energy  expression.  That  is  to  say,  the  last  term  in 
Equation  (2.21)  should  be  retained  to  establish  the  membrane  strain 
energy.  Thus,  if  one  expresses  ( O'/,  (T^' ; ) by  Equation  (2.7) 

with  the  temperature  terms  missing,  the  membrane  strain  energy 

expression  after  the  buckling  will  be  obtained  through  the  following 
procedure 

*5  — ^ £$i'20x<p£x<t>)£>d'>(d4>  4 

'thSff  ( 6j  +2^  e^)&  did<p  -f- 

5 

D Jf  [(  + ^ <~4>)  S £*  t (€<p  £6 4>  + 2( i- o)  £ x<p  S Bdxdtp  + 

5 

1-  D'Jf  C (extDe^)S^  + ie+i  ii6x)Se4+2(i-v)6x4Se»<i>]  Bdijip 
s 

= h Sff  C (7 X 6*  + 0*4  tl(r/+  eU)dd*d<P  -t 

•S 

+ h Sff  ( <r;  ej;  f <%  e l + 2 <rx‘t  e ^ ) 6 d * a * + 

5 

^ 2f  ^ ff  d ^<p  + 2 i>  + 2(1-0)  £>  d zd<p 

u*  = hjf  ( <4  4 + ctq  4 +yr^  e^)5  d*d<$>  + 

+ hff(V*  1-0$  6$  + ZO'x<f>  £x $)%>&'*  d<P  + 
s 

tf  jflel+4f2p€x£*  +2CI-b)6f+]&dvcl4> 

^ 5 


(2.24) 
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The  bending  strain  energy  expression  will  still  remain  Equation  (2.23). 

The  displacement  components  U.",  M",  w"  are  again  assumed,  and 
their  amplitudes  are  also  functions  of  time.  Substituting  the  assumed 
U.''(  V",  in  Equations  (2.6),  (2.8),  (2.19)  and  (2.20),  then  inte- 
grating Equations  (2.15),  (2.23)  and  (2.24)  over  the  total  middle 
• surface  of  the  shell,  the  strain  energy  after  buckling,  TJ=Urn't^b» 
and  the  kinetic  energy  after  buckling,  V » will  be  obtained.  Applying 
the  Hamilton  principle  to  the  new  U and  newV  , a second  order  non- 
linear differential  equation  for  the  generalized  coordinate  of  \aj" 
with  variable  coefficient  is  presented.  This  nonlinear  differential 
equation  will  be  used  to  study  the  dynamic  stability  of  the  shells. 


CHAPTER  3 


LINEAR  SOLUTIONS 


3.1  Conical  Shells 

In  this  section,  Equations  (2.2)  and  (2.4)  are  used  to  substitute 
the  scale  factor  8 and  the  curvature  of  the  shell  In  all  equations 
listed  thereafter  in  Chapter  2.  The  procedure  of  solving  these 
equations  has  been  explained  already  in  Section  2.5.  Now,  the  detailed 
analysis  of  this  portion  will  be  given  in  the  following  sub-sections. 

3.1.1  Transformation  of  coordinates 

According  to  Mushtari  and  Sachenkov  [32],  the  simply  supported 
boundary  conditions  of  conical  shells  can  be  satisfied  by  introducing 


2 = Lo3e  (*/*,) 
or 

y = 'x,  e* 


(3.1) 


Thus,  at  the  ends  of  the  cones  9C=9(|  and  correspondingly,  = q 

an<*  ^0(i«®*»  6 'X2/0 (,) . For  convenience,  transformation  is  also 

made  by  defining 


N/=  <t>Sine 


(3.2) 


According  to  Equations  (3.1)  and  (3.2),  the  following  relations  hold 
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fx  ~(’',e>f'(§j) 

=(«1e5)''(|^)  0.3) 

( 1 Sm9)'2(|^-2)=  (V5)'2^) 

|^[(3<Sln9rl|?] 

By  using  this  transformation,  all  the  basic  expressions  in  the  pre- 
ceding chapter  are  transformed  depending  upon  coordinates  3 and  ^ • 
3.1.2  The  prebuckling  solution 

If  one  assumes  the  temperature  load  T =T (+)  , i.e.,  the 
temperature  is  applied  uniformly  over  the  shell  surface,  and  there 
is  no  temperature  gradient  through  the  thickness  of  the  shell,  then 
one  can  reasonably  expect  that  the  deformation  before  buckling  is 
symmetric.  In  this  case,  as  it  has  been  stated  in  the  preceding 
chapter,  all  the  terms  concerning  v and vanish.  The  displacements 
prior  to  buckling  are  assumed  as 

u'=  Qu(t)Sin  5in(^^)  (3.4a) 

w;=Qw(i)e^5  *‘n(IYL) 

(3.4b) 

in  which  }ra  = loge  A=^,Q(Jtfand  Q*  (t)are  the  ampli- 
tudes of  u'  and  W ' . The  additional  term  ( ( i+  V)o<T  ^09(  5,’h  ( 
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Is  introduced  for  the  purpose  of  satisfying  the  equilibrium  condition 
t27T*ih  +an9CTx]<x  = ^=  ( 27T  *2h  ianQ 

The  method  to  establish  the  strain  energy  and  kinetic  energy 
has  already  been  given  in  Section  (2.5).  The  integration  of  these 
energy  expressions  are  carried  out  over  the  limits  ( J = o to 
^ = o to  ^ = 27T  S*'r»0  ).  The  Hamilton  principle  is  applied  by  partially 
differentiating  the  integrand  with  respect  to  the  generalized  coordinates 
Q'u  and  and  setting  the  results  to  equal  zero  separately.  Thus,  it 
gives  one  algebraic  equation  which  relates  Qu.  and  Qw  by 

F.  Qw  + F3fXT  (3.5) 

and  another  differential  equation  of  with  Qu.  coupled  in  it. 

Considering  Equation  (3.5)  in  this  differential  equation,  the  function 
is  determined  by 

a:  = f,«t  (3.6) 


where  f-  is  a function  of  ( Vj  J»,0  ) and  F2  , depend  on  the  para- 
meters ( v,  R,  , &0, 6 ) • (-jO2"  is  a function  of  ( , R, , &o,  0 , h ) . All 
p p^  are  expressed  in  terms  of  (Ojh;  R ,,  0 , ^ o ) in  Appendix  I. 

Physically,  a)  assumes  the  natural  frequency  and  (/^o< T)  behaves  as  a 
thermal-force.  The  initial  conditions  of  are  reduced  from 
Equation  (2.18)  as 


0. w - Q*  - ° 


at 


t 


The  solution  of  will  depend  on  the  temperature  function  T . 

The  first  type  of  temperature  function  under  study  is  sinusoidal. 
It  assumes 

T = To  Cos  ( cO,t)  (3.7) 
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where  T0  is  a constant  and  a),  is  the  frequency  of  the  temperature 
loading.  The  solution  of  Equation  (3.6)  is  so  obvious  that  it  can 
be  obtained  by  simply  applying  the  Laplace  transform.  It  is  written 
as 


(Q»>S  = -P-V--T  CoS(SL,T) 

s u>z-n .f 


(3.8) 


where  fix  = cO.t  and  T = [ t— ; T. V t . The  subscript  S is  intro- 

LfR(2(l-pz)J 

duced  for  the  case  of  sinusoidal  temperature  function.  Using 
Equations  (3.8)  and  (3.5)  in  Equation  (3.4),  the  prebuckling  dis- 
placements due  to  sinusoidal  temperature  gradient  are  then  determined. 
The  second  temperature  function  is  chosen  as 


T = T,  t/to  for  0 < t < t0 


(3.9) 


T = T,  for  t 

where  T,  is  a constant. 

In  accordance  with  the  conditions  given  in  Equations  (2.18)  and 
(3.9),  the  solution  of  Qw  from  Equation  (3.6)  can  be  obtained  as 
follows.  Let 


X ^ fRfU-V*)} 

Applying  the  Laplace  transform  to  Equation  (3.6)  yields 

Qw  = - S.-nfUt)] 

dQif  _ ot  T,  t Fz  r , _ c oSfcOz)l 
dz  c o2t0  L J 


at  t = t0  or  X - ~ 

t 
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aw  = 

4a'» 

d z 


(S  T,  F2  r t / cOto'i  ] 

to5,n^  JJ 


CKT.tF; 

60 


(3.10a) 

(3.10b) 


When  t >t0  , taking  t - t»  as  the  initial  stage  by  considering 
Equation  (3.10)  as  the  initial  condition,  the  solution  of  (Q,v)r 
is  obtained  by  superimposing  the  particular  solution,  which  again 
can  be  obtained  through  the  Laplace  transform  and  the  complementary 
solution  of  Equation  (3.6),  where  the  subscript  r refers  the 
quantities  to  the  case  of  ramp-type  thermal  loading.  The  explicit 
expression  of  (Qw)r  can  be  written  as 


U)z  v U) 

In  Equation  (3.11) 

Z - 0)(  t - to/t  ) + A 4. 

where  A4  18  a phase  angle. 

For  a Heavyside  load,  t0~°  and 


(3.11) 


(3.12) 


where  the  subscript  h refers  the  quantity  to  the  case  of  a Heavyside 
load. 

Using  Equations  (3.5)  and  (3.11)  or  (3.12)  in  Equation  (3.4), 
the  prebuckling  displacements  due  to  a ramp-type  temperature  loading 
or  a Heavyside  temperature  loading  are  thus  determined. 

The  thermal  stresses  and  moments  are  obtained  by  using  the  above 
solutions  to  the  prebuckling  displacements  in  Equation  (2.7).  Of  course. 
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the  quadratic  terms  related  to  the  normal  displacement,  and  other 
terms  related  to  V and  ^ in  Equation  (2.7)  are  dropped  in  this 
case. 

3.1.3  Considerations  of  dynamic  stability  during  buckling 
During  buckling,  the  additional  displacement  components 
are  chosen  in  the  following  forms 


u"=  Qu  Ct)5.n  ( ggnfcjco 5 (r»,^) 

(3.13a) 

(3.13b) 

w"=  QwCt)  S.o(J!y^)Cos(>l,^) 

(3.13c) 

where  m is  an  odd  integer  and  n,=  n ] SlnO,  in  which  m is  an  integer. 
The  Lagrangian  function  used  in  the  Hamilton  principle  is,  as  it  has 
been  stated  in  Section  2.5,  the  sum  of  Equations  (2.15)  (with 
negative  sign),  (2.22)  and  (2.23).  The  direct  variational  operation 
to  Equation  (2.16)  gives  three  coupled  equations  of  the  generalized 
coordinates  Q*  , and  . Two  of  them  are  algebraic  equations 

If  ))  // 

and  can  be  arranged  to  express  6)  u , Qv  in  terms  of  Ww  . The  third 
one  is  a differential  equation  with  variable  coefficient.  By  con- 
sidering the  first  two  equations,  this  differential  equation  can  be 

_ a 

arranged  to  become  a Mathieu  type  with  only  Qw  as  its  argument. 

For  sinusoidal  temperature  load,  they  are 


(Qu)5=  F+(C&)  $ 

( Os  - 


d*  Os 


Go5f)((Os  = 0 

Ur 


(3.14a) 

(3.14b) 


(3.14c) 
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where 

T = I l,T 


.fl^n  ~ frequency  parameter 

In  general,  F4  and  F^  are  functions  of  ( dD/6,  w,  n ), 

2 

the  frequency  function  _Q„,ndepends  on  ( , Q,  R , /h  , m , n ) and 

the  function  Ttr  is  a function  of  ( }0/  d,  R,/h,  »,  >?/*)  • When 
XL,  = o , i.e.,  oi,  = o , the  smallest  value  of  Ter  with  respect  to 
yn,  r)  assumes  the  static-buckling  temperature  and  is  denoted  by 
(Tcr)Si  • 

For  ramp-type  thermal  loading,  the  final  equations  take 
the  following  forms 


(Qw  — F 4 ( Qw)r 
(ay)r=  F\(Qw)r 


(3.15a) 

(3.15b) 


^■r+^('-^-r-^C0S?)fQwV=O  (3.15c, 

where  Tr  is  a function  of  ( 3.,0,  R,/h,  w,  n,  «.  )•  All  F4  , Fr, 

Tcr  and  Tr  are  expressed  in  Appendix  I.  The  region  of  stability 
is  determined  from  the  parameters  in  Equations  (3.14c)  or  (3.15c). 


3.2  Cylindrical  Sheila 

In  this  section,  Equations  (2.3)  and  (2.5)  are  used  to  substitute 
the  scale  factor  B and  the  curvature  of  the  shell  in  all  equations 
listed  thereafter  in  Chapter  2.  The  procedure  of  obtaining  the  final 
equations  which  indicate  the  characteristics  of  dynamic  stability  of 
this  kind  of  shell  is  parallel  to  that  stated  in  Section  3.1. 
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Therefore,  in  the  following  sub-sections,  results  will  be  simply  given 
without  further  explanation. 

3.2.1  The  prebuckling  solution 

For  a cylindrical  shell,  the  displacement  components  prior 
to  buckling  are  assumed  as 


He  = (Sue  it;  (5in(^)-  if  ) 

Wc  = Qwc  CO 


(3.16a) 

(3.16b) 


where  the  subscript  c.  ifl  used  to  indicate  the  case  of  cylindrical 
shells. 

The  integrations  of  energy  expressions  are  carried  out  over 
the  limit  ( * = ± L/2.  , <t>=o  to  <*>  =27T).  After  the  Hamilton  principle 

is  applied  in  a direct  sense,  the  time-dependent  amplitudes  of 
and  under  a sinusoidal  temperature  loading  are 

(Qnh  = ^ (Q»c)s  (3.17a) 

(0.  U)s  = FytfTo  Cos  (Xl,X)/(tO|  - fit)  (3.17b) 


while  under  a ramp-type  temperature  load  are 

(Qnc)j--  -jj ( ^wc)r  (3.18a) 

( Qwc)r  = £l  aTi  [,.i-  f 2-2  Cos  ) CoS(OcC)  (3.18b) 

assumes  the  prior  buckling  natural  frequency  function.  Both  u)^ 
and  are  functions  of  ( R , h , L , v ) , and  are  expressed  in  Appendix  II. 
For  a Heavyside  load,  to  = o and 
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( ~ Jfk  ( (3.19a) 

( Qm)h  - — -■ ' ( i - Cost  ) (3.19b) 

cot 

The  thermal  stresses  and  moments  can  be  found  in  exactly  the 
same  way  as  in  Section  3.1.2. 

3.2.2  Considerations  of  dynamic  stability  during  buckling 
During  buckling,  the  additional  displacements  are  chosen  in 
the  following  forms 


Vc  = q"c (V  Cos(Z^-)cos(r><P) 

V'Jc.  — Q* C(*J  C0s[^~)  Sin  ( rup) 


(3.20a) 

(3.20b) 

(3.20c) 


where  m is  an  odd  integer,  and  n is  an  integer.  After  applying  the 
Hamilton  principle  in  a direct  sense,  the  final  governing  relationship 
also  assumes  the  type  of  Mathieu's  differential  equation.  For  a 
sinusoidal  temperature  load 


,2  //  . 

d ( Q#c)s 

cItz 


Af  ‘ 


S-  COS  x)( t)s  - O 
•cr 


For  a ramp-type  temperature  load 


(3.21) 


(3.22) 


where  the  superscript  c is  used  to  indicate  the  case  of  a cylindrical 
shell.  In  the  above  equations,  T = gOc(T~  4?  )+A,  A^is  the  phase 
angle,  fl£M  is  a function  of  ( C,  -j-  , m,n  ) and  T<fr  is  a function  of 
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R c 

,n,T1  ).  When  odc equals  zero,  the  smallest  value  of  TCr  with 

respect  to  rn,n  represents  the  static-buckling  temperature  of  the  shell. 

The  function  Tr  also  depends  upon  (-^j-,L,°(,  vn,  n ),  All 

are  expressed  in  Appendix  II.  The  region  of  stability  is  determined 

from  the  parameters  in  Equations  (3.21)  or  (3.22). 

3.3  Numerical  Results  and  Discussions 

The  condition  of  dynamic  stability  depends  upon  the  numerical 
values  of  the  parameters  in  Equations  (3.14c)  or  (3.15c)  for  conical 

shells  and  in  Equations  (3.21)  or  (3.22)  for  cylindrical  shells.  The 

R R 

geometric  parameters  used  for  the  case  of  conical  shells  are  ( h ’ H ' 9 > 
where 

R=(Ri  + R2)/2  0.23) 

H = - 7f()  CoS  0 (3*24) 

The  geometric  parameters  used  for  the  cylindrical  shells  are  R/h  and 
R/L,  where  R and  L are  the  radius  and  length,  respectively.  In  the 
numerical  calculation,  Poisson's  ratio,  , is  taken  as  0.3.  The  length 
of  ramp-type  loading,  t^,  is  assumed  to  b e^li.  For  Heavyside  loading, 
to  is  zero. 

- 200,  500,  1,000;  -g - 0.5,  1.0;  6 - 3#,  10°,  20*; 

200,  500,  1,000;  0.2,  0.5,  1.0 

Values  of  fl,  are  given  as  shown  in  the  figures. 

The  calculations  are  made  by  using  IBM  709.  A flow-diagram  of 
the  Fortran  programme  is  given  in  Appendix  III.  The  parameters 
are  plotted  in  the  diagrams  of  Mathleu's  equation.  In  Figures  2-6, 
the  shaded  areas  represent  the  stable  region.  These  effects 
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of  the  temperature  frequency  and  geometry  of  the  shell  on  the  stability 
are  discussed  below. 

3.3.1  Influence  of  the  force  frequency 

From  Figures  2 and  3,  one  can  notice  that  the  shells  are  more 
unstable  when  the  force  frequencies  -0-,  are  close  to  the  prior  buckling 
natural  frequencies  u)c  and  cO  . This  conclusion  is  similar  to  the 
conclusion  of  Bieniek,  Fan  and  Lackman  [11].  In  reference  [11],  a 
cylindrical  shell  is  considered  under  a harmonic  external  pressure. 
Shells  with  XL,  far  away  from  and  smaller  than  cOc  and  cO  have  not  any 
combination  of  m and  v falling  into  the  first  unstable  region.  This 
can  be  interpreted  that,  for  the  low  force  frequency,  the  stability 
of  a dynamic  problem  is  approximately  the  same  as  the  stability  of  a 
corresponding  static  problem  of  conical  and  cylindrical  shells.  Shells 
with  -XL,  far  away  from  and  greater  than  tOc  and  X)  will  have  some  combi- 
nations of  yn  and  vi  falling  into  the  first  unstable  region.  This  means 
that,  for  the  high  force  frequency,  the  stability  of  a dynamic  problem 
is  quite  different  from  that  of  a corresponding  static  problem  of 
conical  and  cylindrical  shells.  In  fact,  it  tells  that  for  this  case 
the  stability  of  a dynamic  problem  is  much  more  serious  than  the  static 
one.  The  dangerous  phenomenon  of  A(  close  to  tOc  and  u>  can  be  explained 
by  looking  into  Equations  (3.8)  and  (3.17b).  For  a limiting  case,  say, 
A,  = cOc  or  (X)  , then  (Qw)sand  ( Q me  )$  approach  to  infinity.  This  is 
the  case  of  resonance  of  a force  vibration.  Theoretically  speaking, 
the  shells  under  these  conditions  can  never  be  stable  no  matter  how 
small  the  applied  load  is.  Therefore,  one  can  expect  that  force  fre- 
quencies, closer  to  the  resonance  force  frequency,  provide  more  un- 
favorable conditions  to  the  stability  of  the  shell. 
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3.3.2  Influence  of  the  radius-length  parameter 

According  to  Figures  4 and  5,  the  radius-length  parameters 
JL  and  -5.  are  not  important  factors  in  dynamic  stability  problems 
of  cylindrical  and  conical  shells.  This  means  that  for  a certain 
value  of  radius,  thickness  and  apex  angle  (for  conical  shells), 
shells  with  different  lengths  will  become  unstable  at  approximately 
the  same  value  of  temperature  under  a harmonic  thermal  load.  This 
statement  is  more  exact  for  a cylindrical  shell  rather  than  a conical 
shell,  because  as  one  can  see  from  the  figures,  the  effect  of  length 
to  conical  shells  is  larger  than  the  same  effect  to  cylindrical  shells. 
However,  it  seems  that  there  is  a certain  characteristic  difference 
between  conical  and  cylindrical  shells,  since  Figure  5 shows  that  the 
shorter  the  shell  is,  the  easier  it  becomes  unstable,  while  Figure  4 
shows  just  the  opposite.  This  difference  is  physically  interpreted 
by  considering  the  effect  of  the  radius  R,  at  the  upper  end  of  a conical 
shell.  For  the  same  R , a longer  shell  has  a smaller  R,.  This  makes 
it  easier  to  become  unstable.  To  verify  this  intuitive  argument. 

Figure  6 is  plotted  under  the  parameter  . For  the  same  R,,  Figure  6 
shows  that  the  shorter  shell  becomes  unstable  first.  This  agrees  with 
the  trend  of  the  cylindrical  shell.  For  the  ramp-type  temperature  load, 
the  effect  of  the  radius-length  parameter  is  also  very  small. 

3.3.3  Influence  of  the  radius-thickness  parameter 

From  Figures  4,  5 and  6,  one  can  immediately  observe  that  the 

R R 

radius-thickness  parameters  and  are  the  most  important  geometric 
factors  for  dynamic  stability  problems  of  conical  and  cylindrical  shells. 
These  figures  show  that  for  the  same  value  of  the  radius-length  para- 
meter and  the  apex  angle  (for  conical  shells),  different  values  of 
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JL  or  Ji  will  influence  the  stability  of  the  shell  to  different  extents, 
h h 

a f? 

The  larger  the  value  of  -7-  or  -r  is,  or  say,  the  thinner  the  9hell  is, 

n h 

the  more  unfavorable  the  condition  of  dynamic  stability  of  the  shell 
will  be.  This  conclusion  is  in  accordance  with  the  conclusion  of  the 
static  case.  For  a ramp-type  temperature  load  or  Heavyside  temperature 
load,  the  conclusion  is  also  the  same  as  above.  A series  of  curves  are 
plotted  in  Figures  7 and  8. 

3. 3. A Influence  of  the  apex  angle 

R o 

In  Figure  A,  for  each  set  of  values  of  -1-  and  , three  differ- 

H n 

ent  values  of  the  half-apex  angle  6 are  imposed.  The  curves  show  that 
the  apex  angle,  just  as  the  length,  plays  an  unimportant  role  in  a 
dynamic  stability  problem  of  a conical  shell.  However,  it  reveals  that 
the  shell  under  a harmonic  temperature  load  will  be  more  unstable  as  the 
apex  angle  is  increased.  For  a ramp-type  temperature  load,  the  con- 
clusion is  again  the  same  as  above. 

3.3.5  Comparison  between  the  effect  of  a ramp-type  temperature 
load  and  a Heavyside  temperature  load 

As  it  has  been  mentioned  before,  the  Heavyside  temperature 
load  can  be  deduced  from  the  ramp-type  load  simply  by  setting  t0=°. 
Therefore,  in  Figures  7 and  8,  curves  with  = o and  te  = Tit  are 
plotted  separately  with  T,  equal  to  the  static-buckling  temperature 
of  the  shell  under  uniform  temperature  load.  The  results  show  that 
the  Heavyside  temperature  load  i9  more  unfavorable  than  the  ramp-type 
temperature  load.  These  results  are  also  similar  to  the  results  of 
a cylindrical  shell  under  the  same  types  of  external  pressure  [11]. 

In  fact,  these  results  can  be  expected  intuitively.  Since,  if  one 
transforms  the  origin  of  the  time  coordinate  from  t ~o  to  t =ts  for 
the  ramp-type  load,  one  will  find  that  shells  under  a ramp-type  load  are 
actually  the  same  as  shells  under  a Heavyside  load  with  initial  de- 
flection and  velocity. 


CONVERSION  OF  COEFFICIENTS  OF  THE  LINEAR  MATHIEU  DIFFERENTIAL  EQUATION 
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Figure  2.  Mathieu's  Diagram  of  Conical  Shell  under  Sinusoidal 
temperature  Loading  with  Three  Frequencies  _Q.,=o.3  , i.c 

\-5  ( R/H  = l , cxT0  = 3.5"  X 1 0 , cO  = o-  77  , 6=20°  ). 
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Figure  3.  Mathieu's  Diagram  of  Cylindrical  Shell  under  Sinusoidal 
Temperature  Loading  with  Three  Frequencies  .Q. , = o • 3 , 

io  ,j-5<  (R/l=I.  cXTo  = °<3  X /o'+  , Gdc=o.?<?). 
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Figure  4.  Mathieu's  Diagram  of  a Conical  Shell  under  Sinusoidal 
Temperature  Loading  ( XL,  = I ■ !T  , o<To  = 3. 5 x |0  "4  ). 
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Figure  5.  Mathieu's  Diagram  of  a Cylindrical  Shell  under 
Sinusoidal  Temperature  Loading  , «.T0  = 3-5  X 10*4  ). 
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Figure  6.  Mathieu's  Diagram  of  a Conical  Shell  under  Sinusoidal 
Temperature  Loading  ( jt,  = 1.5  , a T„  = 3-5  A I0~4  , 9-20°  ). 


rxio 

Figure  7.  Mathieu's  Diagram  of  a Conical  Shell  under  Ramp— Type  and  Heavyside 
Temperature  Loadings  ( , 9 = 2 0°  , T./iTcr)*  = I ). 
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Figure  8.  Mathieu's  Diagram  of  a Cylindrical  Shell  under  Ramp-Type 
and  Heavyside  Temperature  Loadings  ( L = 2R  , T,  /(TciOs-t  = I )• 


CHAPTER  4 


% 


NONLINEAR  SOLUTIONS 

4.1  Conical  Shells 

In  this  section,  just  as  in  Section  3.1,  the  scale  factor  B 
and  the  curvature  of  the  shell  1?^,  appearing  in  all  basic  equations 
of  Chapter  2,  are  expressed  by  Equations  (2.2)  and  (2.4).  The  co- 
ordinates are  transformed  to  the  pair  of  coordinates  ( Jf  ^) 

through  the  transformation  explained  in  Section  3.1.1.  The  detailed 
analysis  of  the  conical  shell  through  the  procedure  stated  in 
Section  2.5  is  given  in  the  following  sections. 

4*1*1  Basic  governing  nonlinear  differential  equations 

The  prebuckling  solution  to  the  present  case  is  exactly  the 
same  as  obtained  in  Section  3.1.2.  The  prebuckling  deformation  is 
assumed  to  be  symmetric.  After  buckling,  the  additional  displace- 
ment components  are  chosen  in  the  following  forms 


U"  = A,(  +)S,'n  ( —^)CoS(n^)-^  A2(t)5m(A^;c0S(2^»  + 

1-  AjltjS.n  f (4.1a) 

d o ' 


v"=  B,Ct)e 

(— y ) S in  ( + 

+ 02^)  5 Wi 

(2n,*)  [CoS(2n^*)  l) 

(4.1b) 

\aj  "~  c,u;  e^sin 

(4.1c) 
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where  A,(t),  A2(t),  A^*);  B,(+>,  B-jlOand  (^(+)  are  time-dependent  ampli- 
tudes of  u" , \/"and  w".  The  additional  terms  ( ^HlJL - JCoffityf') 

A3  (t)  5 in  ( - an<*  6 ^j^)-jJs/ri(2'iiyare  introduced 

for  the  purpose  of  showing  the  effects  of  the  nonlinear  coupling  terms 
in  the  strain  energy  expression.  If  Equation  (2.24)  is  expressed  in 
terms  of  the  displacement  components  u",  \t"  and  tv",  it  will  assume 
the  following  form  for  the  conical  shell. 


u = 


r3o  r2VS;»d 

h II  ( (Tx  U.)  + Of  (y,+  - c oi  0w"f  U"J 


o o 


+ Or* $ ( v,}  -t  u%  -v  ")](*,  e*)d$d  + + 


. r%o  r JVSir,  9 

+ 2 J0  jo  C °}7 v*  ( *,t)z+  2 crx;  ( ^)(  w/jjj  fydij,  -f 

, rZt,2VSm9 

+ 2j  J ((  *"$)Z+ ( - C0tGw"-t  u")Z+  ZvU.^  ( \l’y  - Coie  tv'V 


+ u.")  + ~ r + u%  - o27 ^ ^ + 


^ i to @ 

+ J]  [(*%)*  + (Wfy4  + 2(w, If (W,”+)1]  (2*le*)~2d}d'l' 

O J o 


+i 


/ rhf2^9 

j ]o  f usZ(w,l)z  + ( v,+  - cole  a")(*j;f)2+  u u^(^l)z-f 


t ^ f - co+d  w"-t  u")( tv"/  + ( / -m) r 1/ " t u'j^, - u "5 tv'' iv/^J 


(* .el)-' d & dt 


(4.2) 
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If  the  additional  terms  are  dropped  from  Equations  (4.1a)  and  (4.1b), 
one  can  immediately  observe  that  after  Equations  (4.1a),  (4.1b)  and 
(4.1c)  are  substituted  into  Equation  (4.2),  the  last  integral  of 
Equation  (4.2)  equals  zero  after  integration.  If  this  is  the  case, 
then  there  would  not  be  any  relaxation  to  the  nonlinear  spring  force 
which  is  introduced  by  the  fourth  integral  in  Equation  (4.2).  In 
other  words,  the  tangential  displacement  components  u."  and  v*  do  not 
play  any  role  in  restraining  the  enlargement  of  the  normal  displacement 

II 

component  W , even  the  ends  are  completely  restricted.  Apparently, 
this  is  not  true.  Therefore,  the  additional  terms  of  Equations  (4.1a) 
and  (4.1b)  must  be  introduced  at  this  stage. 

The  Lagrangian  function  used  in  the  Hamilton  principle  is, 
as  it  has  been  described  in  Section  2.5,  the  sum  of  expressions  of 
Equations  (2.15)  (with  negative  sign),  (2.23)  and  (2.24).  The  direct 
variational  operation  to  Equation  (2.16)  gives  six  equations  of  the 
generalized  coordinates  A,,  A1;  A?  . 0,  , and£(  . Because  of  the 
orthogonality  among  certain  terms  of  Equation  (4.2),  A,  and  ^ are 
uncoupled  with  and  62  , while  is  uncoupled  with  all  A,  , Ai;g| 
and  B2 . Five  of  the  six  equations  are  algebraic  ones;  thus,  they 
can  be  solved  simultaneously  to  express  ft1(  and  B^in  terms 

of  ♦ The  las t one  is  a nonlinear  differential  equation  with  vari- 
able coefficient.  By  considering  the  first  five  equations,  this 
differential  equation  can  be  arranged  to  have  its  only  argument. 

« 
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For  sinusoidal  temperature  load,  they  are 


<fl,ls=F  4(d,)f 


(4.3) 


(6,),=  F5(C,)s 


(4.4) 


(4.5) 


(b2)s=fs(5,)2(-£;) 


(4.6) 


in  which  the  bar  over  the  generalized  coordinates  F\hF\Zl  A3  ; B,,  8Z 
and  C,  indicates  that  these  generalized  coordinates  are  nondimension- 


yl(,yl2.&nd  are  functions  of  ( 3o,9,>n,  E ) . \ is  the  nonlinear 

force  introduced  by  the  first  term  of  U.1'  and  \j"  . is  the  nonlinear 
force  introduced  by  the  second  term  of  H."and  y"  , while  ^3is  the 
nonlinear  force  introduced  by  the  third  term  of  u" . They  are  listed 
in  Appendix  I . 

For  ramp-type  thermal  loading,  the  final  equations  take  the 
following  forms 


In  general,  Fg , Fg  and  F(0  are  functions  of  ( £<,,0,  yn,  n ) , and 
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(Mr  = P4(C.)r 


(4.9) 


(6.)r  = F5(C',)r 


(4.10) 


(4.13) 


(4.12) 


(4.11) 


It  is  appropriate  to  point  out  here  that  if  all  nonlinear 
terms  of  (C()5  and  ( C,  )r  are  dropped,  then  the  above  two  sets  of 
equations  reduce  to  the  corresponding  linear  ones.  As  one  may 
observe,  they  are  exactly  the  same  as  those  obtained  in  the  pre- 
ceding chapter  [Equations  (3.14a),  (3.14b)  (3.14c);  (3.15a),  (3.15b) 
and  (3.15c)]. 

Equations  (4.8)  and  (4.14)  are  the  governing  equations  for 
studying  the  problem  of  nonlinear  dynamic  behavior. 

4.1.2  Nonlinear  vibrations 

If  terms  which  include  temperature  T»  or  T,  in  Equations  (4.8) 
and  (4.14)  are  set  to  zero,  i.e.,  if  the  shell  is  free  of  any  external 
temperature,  then  both  equations  reduce  to  the  following  unique  form. 


(4.15) 
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Equation  (4.15)  governs  the  nonlinear  free  vibration  of  a conical 
shell.  The  solution  of  Equation  (4.15)  can  be  expressed  as 


oo 

C,  = Q;  ( COS  Cl' SIT  )+  5m(«'iU  )) 

l-l 


(4.16) 


As  a first  approximation,  £,  can  be  written  as 


C,  - P,(  CoS  (SIX)  + Si'n(Slz))  ( 

Substituting  Equation  (4.17)  into  Equation  (4.15)  and  comparing 
coefficients  of  corresponding  modes  after  calculation  yields 


(4-18> 

Equation  (4.18)  gives  the  relation  between  the  frequency  of  free 
vibration  n.  and  the  corresponding  vibrating  amplitude  a,. 

4.1.3  Dynamic  stability  and  response 

One  may  notice  that  both  Equations  (4.8)  and  (4.14)  are 
of  the  same  form  as  the  so-called  nonlinear  Mathieu  differential 
equation  [33].  For  the  purpose  of  simplicity,  the  following 
Equation  (4.19)  is  used  to  develop  the  solution.  One  Just  compares 
Equations  (4.8)  and  (4.14)  with  Equation  (4.19);  then  the  corresponding 
solutions  will  be  obtained. 

jfz  =0  (4.19) 

It  has  been  known,  [2],  [30],  [33]  and  [34],  that  the  in- 
stability regions  of  the  linear  solution  of  Equation  (4.19)  are 
located  near 


-CL  ~ 2.  -TL  inn  / £ 


( 4 = ',  2,  - ) 


(4.20) 
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The  most  dangerous  of  these  regions  is  the  first  instability 
region  which  has  been  used  in  discussing  the  dynamic  stability  for 
the  linear  case  in  Chapter  3.  Therefore,  as  a first  approximation, 
the  solution  of  Equation  (4.19)  will  be  taken  in  the  following  forms 
which  are  near  the  principal  (first)  instability  region  [30]. 


C,  = b,  COS  1 
/ / 2 


(4.21a) 


for  the  branch  adjoining  the  lower  boundary  of  the  instability  region, 
and 


Ct  - C,  Si n jt 
/ 2 


(4.21b) 


for  the  branch  adjoining  the  upper  boundary  of  the  instability  region. 

In  Equation  (4.21),  b(  and  c(  are  two  constant  coefficients  representing 
the  amplitudes  of  the  normal  displacement  component  of  the  shell.  After 
substituting  Equations  (4.21a)  and  (4.21b)  into  Equation  (4.19),  b,  and 

C,  can  be  solved  through  the  process  of  comparing  coefficients.  They 
are 


± (i_2 

3/3  1 4 


+M  ) 


(4.22a) 


l-M  ) 


(4.22b) 


where  ^2-  These  expressions  coincide  with  the  corresponding 

relationships  for  plates  [2]. 


It  is  of  interest  to  point  out  that  as  long  as  one  remembers 
that  the  term  containing  Cos 't.  or  Sfn  is  the  dominant  term  in  the 
solution,  a second  approximation  of  the  solution  gives  exactly  the 
same  results  for  b(  and  c,  as  obtained  by  the  first  approximation. 

As  a second  approximation,  the  solution  of  C,  will  be  taken  in  the 
following  forms 
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C'/=  bo  + t>,C0S±  +b2C0S$  (4.23a) 

c,  = Ce  -h  c,Sm  | 4-  C2  S//10-  (4.23b) 

The  reason  for  inclusion  of  the  free  terms  bc  and  c,  is  self-evident. 
Indeed,  inserting  Cos  y into  the  nonlinear  part  of  Equation  (4.19), 
one  will  arrive  at  an  expression  with  a constant  term  of  the  same 
order  of  magnitude  as  the  coefficient  of  Cos Inserting  Equation  (4.23a) 
into  Equation  (4.19),  and  taking  b,  » ba  , b,  » b2  , one  will  obtain 
three  nonlinear  algebraic  equations  containing  be/ b,  and  bz  by  equaling 


the  coefficients  of  the  free  term,  CoS  1 and  CoSJ  to  zero, 
equations  are  of  the  following  forms 

These  three 

(1  ■+  f-  p bf ) ba  -t  ( ^ fib?  ) b2  ss  0 

(4.24a) 

( j/3b/2-2/uj b0  + ( J-zSb2 f / - t,2)  b2  = ° 

(4.24b) 

(4.24c) 

Equation  (4.24c),  as  one  can  notice,  is  exactly  the  same  as 
Equation  (4.22a).  Thus,  Equations  (4.24a)  and  (4.24b)  give 
only  a trivial  solution. 

b0  and  bz 

the  same 

Of  Co,  C 

Inserting  Equation  (4.23b)  into  Equation  (4.19)  and  making 
simplification  c,»c0  » C,  » Cz  . three  nonlinear  equations 
n cz  will  be  obtained  as 

( 1+  %/6C?)  C 0 - ( M + ^/3C,Z)  C2  = o 

(4.25a) 

~(*M+\lbC*)C0  -Ki-£,2+  f/3 cz,)Cz~o 

(4.25b) 

\ + 'u  + 4 /^N0 

(4.25c) 
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Equation  (4.15c)  is  the  same  as  Equation  (4.22b).  Thus,  Equations 
(4.25a)  and  (4.25b)  give  C0  and  only  a trivial  solution. 

From  Equations  (4.24c)  and  (4.25c),  one  may  plot  bf  and  C, 
against  the  frequency  parameter  and  thus  form  the  so-called 
response  curves  in  the  amplitude-frequency  plane.  Usually,  the 
absolute  value  of  b(  and  C,  is  used  as  the  ordinate,  while  is 
used  as  the  abscissa  (See  Figs.  9a,b,c).  When  /2>~o>  these  response 
curves  reduce  to  their  counterparts  of  a linear  forced  vibration 
(See  Fig.  9c).  When  /(  = o,  Equations  (4.24c)  and  (4.25c)  become 
identical.  The  response  curve  for  M-°  is  drawn  heavier.  Com- 
paring Figures  9a, b with  Figure  9c,  one  sees  that  the  response  curves 
in  the  nonlinear  cases  could  be  thought  of  as  arising  from  those  for  4 

the  linear  case  by  bending  the  latter  to  the  right  for  a hard  spring 
and  to  the  left  for  a soft  spring.  Also  from  these  two  figures,  one 
sees  that  for  the  linear  vibration  the  relation  between  the  ampli- 
tude and  the  frequency  is  conformal,  while  for  the  nonlinear  case 

A 

it  is  not.  For  certain  values  of  ^ , there  could  be  two  different 
values  of  b or  C(  . 

Because  of  the  nonconformal  relation  between  the  amplitude 
and  the  frequency,  the  so-called  jump  phenomenon  happens.  This 
phenomenon  will  provide  the  information  of  the  stability  condition  of 
the  shell  under  consideration.  Let  us  imagine  an  experiment  performed 
in  which  M is  held  constant,  while  si  is  slowly  varied,  and  the  ampli- 
tudes b and  c,  of  the  resulting  harmonic  vibration  are  observed. 

Consider  first  the  case  of  a hard  spring  force,  /3>o,  and  suppose  that 
4"  is  rather  large  at  the  beginning  of  the  experiment,  i.e.,  by  re- 
ferring to  Figure  10a,  it  is  started  at  point  1 on  the  curve.  As  ^2is 
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&>0  (Hard  Spring)  fi  < 0 (Soft  Spring) 


(3) 


(b) 


& -o 


(Linear  Case) 

(C) 


Figure  9.  Amplitude-Frequency  Diagram 
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decreased,  C,  immediately  jumps  to  fc>(  at  point  2,  then  follows  the  path 
through  point  3.  Upon  performing  the  experiment  in  the  other  direction, 
i.e.,  starting  at  point  4 and  increasing  ^2,  the  amplitude  follows  the 
path  4-3-2  to  infinity  without  any  jump.  In  the  real  case,  there  are 
always  damping  forces  in  the  shell  no  matter  how  small  they  are.  There- 
fore, for  the  present  case,  one  may  imagine  that  the  solution  of  the  ampli- 
tude following  the  path  4-3-2  will  jump  at  a point  5 near  infinity  to 
point  6 on  the  abscissa.  The  circumstances  are  quite  similar  with  a 
soft  spring  force,  but  the  jumps  in  amplitude  take  place  in  the  reverse 
direction  (See  Fig.  10b). 


/3  > O 

(a) 


Figure  10.  Jump  Diagram 


It  seems  plausible  that  the  solution  corresponding  to  absolute  values 
of  the  amplitude  and  the  frequency  parameter  for  those  portions  of  the 
harmonic  response  curves  which  are  passed  over  in  jumping  from  one 
point  to  another  are  unstable  solutions.  In  other  words,  for  the  hard 
spring  force,  the  solution  originating  at  the  upper  boundary  of  the 
instability  region  is  unstable,  while  for  the  soft  spring  force,  the 
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solution  originating  at  the  lower  boundary  of  the  instability  region 
is  unstable.  This  is  shown  in  Figure  10. 

In  discussing  conical  shells,  one  may  compare  Equations  (4.8) 
and  (4.14)  with  Equation  (4.19)  to  find  out  their  correspondences. 
For  the  sinusoidal  temperature  load.  Equations  (4.24c)  and  (4.25c) 
take  the  following  forms 


_l_ 

4 


n: 


n 


mn 


2 Tc°r+  ^ ^3)b f-o 


' mn 


(4.26a) 


4 sii„ 


+ 4_i  + 

2 Ter 


(4.26b) 


For  the  ramp- type  temperature  load,  Equations  (4.24c)  and  (4.25c) 
take  forms  as 


o)Z 

Tcr  4 nL 


2 Tr+4  ^ b'  " 


0 


(4.27a) 


J_  cO2  J_  T, 

4 2 Tr 
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-(■ 
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■O.Z, 


^~Vcf  =° 


(4.27b) 


The  correspondence  between  coefficients  of  Equations  (4.24c)  and  (4.25c) 
and  Equations  (4.26a),  (4.26b);  (4.27a),  (4.27b)  is  shown  in  Table  II. 


CONVERSION  OF  COEFFICIENTS  OF  THE  NONLINEAR  MATHIEU  DIFFERENTIAL  EQUATION 
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4.2  Cylindrical  Shells 

In  this  section,  the  scale  factor  B and  the  curvature  of  the 
shell  1^,  appearing  in  all  basic  equations  of  Chapter  2,  are  expressed 
by  Equations  (2.3)  and  (2.5).  The  detailed  analysis  of  the  cylindrical 
shell  through  the  procedure  stated  in  Section  2.5  is  given  in  the 
following  sections.  Since  the  general  approach  is  the  same  as  that 
for  a conical  shell,  certain  steps  which  have  been  discussed  thoroughly 
in  the  above  section  will  not  be  repeated  in  the  present  section  in 
detail. 

4.2.1  Basic  governing  nonlinear  equations 


as  those  in  Section  3.2.1.  All  equations  in  Section  3.2.1  are  used 
for  further  analysis  in  the  present  section. 


The  prebuckling  solutions  to  the  present  case  are  the  same 


After  buckling,  the  additional  displacement  components  are 


assumed  as 


+ S;n(^j^)  Cos(ZTi<t>)  + 


+ AKCt)  Sfn(^I2L) 


(4.28a)  • 


Vc  = COS  ( 22J*)  CoS  + 


+ B2t(*)  ( Cos(^)  ■+  \ )S!»  (2n4>) 


(4.28b) 


Wc  = Clc(t)  cos 


(4.28c) 
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The  additional  terms  Alclt)  Sin  ) CoS  (zn<|>),  A 3C  (<")  S.'n(— ~) 

and  B?cLt)(  Co$(  + | ) Sin  (2  n 4>)  are  introduced  for  the  same  purpose 

as  explained  in  the  last  section,  i.e.,  they  are  introduced  to  show 
the  effects  of  the  nonlinear  coupling  terms  of  the  displacement 
components  in  the  strain  energy  expression. 

The  Lagrangian  function  used  in  the  Hamilton  principle  is 
the  sum  of  expressions  of  Equations  (2.15)  (with  negative  sign),  (2.23) 
and  (2.24).  The  direct  variational  operation  to  Equation  (2.16)  gives 
six  equations  of  the  general  coordinates  A<c/  A3C  ’»  B,c  , B2cand 
CIL.  Through  the  same  arrangement  described  in  Section  4.1.1,  the 
final  governing  nonlinear  differential  equations,  corresponding  to 
Equations  (4.8)  and  (4.14)  for  the  conical  shell,  assume  the  following 
forms 


dZ(cl(.)s  + (nL)2^_  cos  z)( CICK 


d x‘ 


-a. 


Tc 


cr 


°?ic~  ~^l)c  (j, 

C,t*S  " (4*29) 


d2(C,c)r  +<*L),  , T.  _ JL  o 

co*  (/  Ter  rr^Sx)(C,t)r  co*  ^ (4*30) 


where  and  are  functions  of  ( ^ >*7,  >1  ) (See  Appendix  II). 

If  nonlinear  terms  in  Equations  (4.29)  and  (4.30)  were  dropped,  these 
two  equations  also  reduce  to  Equations  (3.21)  and  (3.22)  obtained  in 
the  preceding  chapter. 

4.2.2  Nonlinear  vibrations 

Similar  to  the  conical  shell,  if  temperatures  To  and  T,  are 
dropped  from  Equations  (4.29)  and  (4.30),  then  these  two  equations 
reduce  to 

77?  + ' ) ■ 2 c,t  ( x-  v - n.  ) c?  = o 


(4.31) 
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Assuming 


C,c  - Q,c  Cos  (Hz) 


(4.32) 


and  substituting  Equation  (4.32)  into  Equation  (4.31)  yields 


Equation  (4.33)  relates  the  frequency  of  the  free  vibration,  TL  , and 
the  corresponding  vibrating  amplitude  9(C . 

4.2.3  Dynamic  stability  and  response 

Equations  (4.29)  and  (4.30)  assume  the  same  form  as 
Equation  (4.19),  Therefore,  they  are  also  governed  by  the  criterion 
developed  in  the  latter  part  of  Section  4.1.3. 

For  the  sinusoidal  temperature  load,  Equations  (4.24c)  and 
(4.25c)  take  the  following  forms 


For  the  ramp-type  temperature  load,  Equations  (4.24c)  and 
(4.25c)  take  forms  as 


(4.33) 
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The  correspondence  between  coefficients  of  Equations  (4.24c),  (4.25c) 
and  Equations  (4.34a),  (4.34b)  and  (4.35a),  (4.35b)  is  also  shown  J.n 
Table  II. 

4 . 3 Numerical  Results  and  Discussions 

For  additional  interest,  certain  numerical  results  for  both 
cylindrical  and  conical  shells  are  calculated  through  IBM  709  and 
are  plotted  into  curves  for  several  governing  parameters.  Values  of 
geometric  parameters  are  chosen  to  be 

-y  - 500,  ~ - 0.5,  9 ■=  10° 
for  a conical  shell,  and 

-5-  - 500,  - 0.5 

for  a cylindrical  shell. 

The  force  frequency  il, , is  chosen  to  equal  2Hmnor  3 -ft-*,,, 
(jfl^for  cylindrical  shells).  Poisson's  ratio,  v , is  taken  as  0.3. 
to  is  assumed  to  be  JTt  for  the  ramp-type  temperature  loading,  and 
zero  for  the  Heavyside  temperature  loading.  The  governing  parameters 
are  picked  up  to  be  il^,n,o<(  Tee  )j  ,o((  Tor  )r  , o(Tr,  ^ and  (A^,„)  ,°C(Ttr)St 
c((Tcr)r,  OCT where  and  ^ zC~  V.  ,c  • 

The  subscripts  s and  r indicate  the  quantities  corresponding  to  sinus- 
oidal, ramp-type  or  Heavyside  temperature  loading.  These  governing 
parameters  are  plotted  against  the  mode  numbers  and  n.  through 
Figures  11  - 20. 

The  free  vibration  of  a conical  shell  for  several  different  modes 
is  shown  in  Figure  21,  while  the  free  vibration  of  a cylindrical  shell 
for  several  different  modes  is  shown  in  Figure  22.  They  are  of  the 
same  form  as  the  free  vibration  of  a simply  supported  plate  [35],  It" 
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is  observed  from  Figures  21  and  22,  that  under  the  same  value  of 

the  larger  the  mode  number  is,  the  smaller  the  vibrating 
amplitude  will  be.  Also  in  Figures  21  and  22,  it  reveals  that  the 
additional  terms  added  to  the  displacement  components  uv,  V"  and 

, V do  influence  the  vibration  to  a certain  degree.  Nevertheless, 
it  seems  that  the  third  term  of  u.  or  Uc  (the  symmetric  term)  does 
a minor  job.  So  it  would  be  of  no  harm  if  it  is  neglected  in  the 
calculation. 

From  Figures  14  and  20,  it  has  been  noticed  that  the  nonlinear 
spring  force  ^ or  is  always  positive.  This  means  that  for  a 
moderate  shell  under  the  present  boundary  conditions,  the  nonlinear 
effect  always  causes  a hard  spring  force.  Therefore,  the  solution 
originating  from  the  upper  boundary  is  always  unstable  while  that 
originating  from  the  lower  boundary  is  always  stable.  This  is  shown 
in  Figures  23  and  24.  Figures  23  and  24  also  reveal  that  //-  and  ft 
are  two  control  parameters  of  the  stability  of  the  shell.  The  width 
of  the  instability  region  bounded  by  Equations  (4.24c)  and  (4.25c)  is 
controlled  by  /A-  , while  the  slope  of  the  boundaries  is  controlled  by 
/3  . Thus,  from  Table  II  combining  with  Figure  23,  one  will  realize 
that  the  higher  the  temperature  gradient  T»  or  T,  is,  the  more 
dangerous  the  stability  condition  for  both  conical  and  cylindrical 
shells  will  be. 

Figure  25,  the  stable  and  unstable  solutions  of  the  conical 
shell  subjected  to  sinusoidal  temperature  loading  with  w=*5, 

~ SO  c aRd  JT  t -3-flrri^  and  those  subjected  to  ramp-type  or  Heavyside 
temperature  loading  with  m = 9 , 71  = 17  and  T(=50°c  have  been  depicted. 

^*,A‘-and£  were  found  by  using  Table  II  and  Figures  11  - 15.  The 
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points  P(  and  pH  are  the  stable  and  unstable  solutions  under  a 
sinusoidal  temperature  loading.  Similarly,  the  points  P2  and  P2i 
are  stable  and  unstable  solutions  under  a Heavyside  temperature  load, 
and  the  points  Pj  and  P33  are  stable  and  unstable  solutions  under  a 
ramp-type  temperature  loading.  By  looking  at  P2  and  P3  , one  comes 
to  the  same  conclusion  that  was  made  for  the  linear  case,  i.e.,  the 
Heavyside  temperature  loading  is  more  unfavorable  than  the  ramp-type 
temperature  loading.  The  stable  and  unstable  solutions  of  the  cylindrical 
shell  can  be  obtained  through  the  same  procedure  by  using  Table  II  and 
Figures  16  - 20.  The  results  are  similar  to  those  of  the  conical  shell. 
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Figure  11.  Frequency  Parameter  Diagram  of  a Conical  Shell 

( R/H  = 0.5,  R/h  = 500  , 0 =10°  ). 
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Sinusoidal  Temperature  Loading  (.R/h=0.5,  R/h -500  k Q-io°  , 

= 4 n mn  , 
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Figure  14.  Temperature  Parameter  Diagram  of  a Conical  Shell 
under  Ramp-Type  and  Heavyside  Temperature  Loadings 

( R/H  = o.5  , R/h  =500  , 0 = 10°  , t»  = Tf  t ,o  ). 


n 
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Conical  Shell  ( r/h  = o.5  , R/h  = 500  , e=io°  ). 
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Figure  16.  Frequency  Parameter  Diagram  of  a Cylindrical  Shell 

( R/l=  0-5  , R/h  = 500  ) 
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Figure  17.  Critical  Temperature  Diagram  of  a Cylindrical  Shell 
under  Sinusoidal  Temperature  Loading  ( R/l=o-5  , 

R / h “ 5 00  , XI ! — 
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Figure  18.  Critical  Temperature  Diagram  of  a Cylindrical  Shell 
under  Ramp-Type  or  Heavyside  Temperature  Loading 

( R/L  = o-S  , R/h  = 500  ). 
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Figure  20.  Nonlinear  Spring  Force  Parameter  Diagram  of  a 
Cylindrical  Shell  ( R/l-o.5  , R /h-500  ). 
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Figure  21.  Nonlinear  Free  Vibration  of  a Conical  Shell 

( r/h=  0 5 , R/h  =500  , 0 = /o°). 
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Figure  23.  Amplitude-Frequency  Diagram  with  Constant  $ . 
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Figure  24.  Amplitude-Frequency  Diagram  with  Constant  ^ . 


r 1.074 


79 


N 


.a- 


Figure  25.  Amplitude-Frequency  Diagram  of  a Conical  Shell 

( R/H=  0.5  , R/h  =500  , 6 -10°,  rtTo^fc-feXIO"^ 

(XT.=  fc.fc  X |0~*  , f.  = 0»7Tt) 


CHAPTER  5 


CONCLUSION 

An  energy  approach,  to  analyze  the  dynamic  behavior  of  conical 
and  cylindrical  shells  under  surface  temperature  load  by  linear  and 
nonlinear  shallow  shell  theories  through  steps  stated  in  Chapter  2, 
is  able  to  find  out  some  characteristics  of  this  kind  of  problem. 

For  the  linear  case,  it  tells  us  that  the  essential  effects  on  dynamic 
stability  of  conical  and  cylindrical  shells  come  from  the  inertia 
terms  of  the  prior  buckling  equation  of  motion  as  well  as  the  radius- 
thickness  ratio  of  the  shell.  Lengths  and  appex  angles  are  less 
influential  factors  in  the  present  problem.  Generally  speaking,  the 
dynamic  stability  is  more  dangerous  than  the  static  stability  under 
the  proposed  boundary  and  loading  conditions.  For  the  nonlinear  case, 
it  reveals  that  additional  terms  of  the  tangential  displacement  com- 
ponents are  needed.  Otherwise,  the  nonlinear  force  calculated  would 
be  larger  than  it  actually  is.  For  shells  with  moderate  size,  under 
the  present  boundary  conditions,  the  nonlinear  effect  always  exercises 
a hard  spring  force.  Thus,  the  response  solution  originating  at  the 
upper  boundary  of  the  instability  region  is  always  unstable,  while 
the  response  solution  originating  at  the  lower  boundary  of  the  in- 
stability region  is  always  stable. 

For  both  conical  and  cylindrical  shells,  the  dynamic  stability 
condition  becomes  more  dangerous  when  the  temperature  gradient  To  or  T, 


80 


81 


is  higher.  The  free  vibration  which  has  infinite  modes  of  frequency 
is  different  from  the  linear  free  vibration  for  which  the  frequency 
of  the  vibration  equals  to  the  natural  frequency  of  the  shell.  For 
sinusoidal  temperature  loading,  higher  temperature  causes  larger 
response.  For  ramp-type  and  Heavyside  temperature  loadings,  the 
latter  causes  larger  response. 
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APPENDICES 


APPENDIX  I 


FORMULATIONS  OF  THE  CONICAL  SHELL 

I = 7T/3-o,  SO  = Sin  0y  K9=Coie 
£/=  exp(2o)-l,  e2  = exP(2\^0)  -/,  e3  = exp  ( \l0)  + i; 

expCr/4A)30J  + i,  E5=  exp  (2^-1)^]-!,  p6=  exp  [(2X-i)^0]-j/ 

E7  = exp  ((3A-i)3<>] -t  \,  Eg  = exp  (r  2 + A;^a]  + 1 ^ E?  - exp  (2  dfX^0J/ 

E, 0=  1 - expt(4K-2)io), 

P —,  + 16 1*  P,=  I2+  A2,  P/(  = I2W,+  A)2  P/((=i2  + (,-A)  2 

P2  = 9i2+A2  P2Z=9lV>+A)z  p222=  9I2+(|_A)2 

A,  = f2I24  1/2)^  A2-  81(1  + 41  )E,P“'  A3  = I2£2/(4AP,); 

A4=  0/P2-  l/P()E3A;  A?  = ( </Pn-  «/P|l)E4n+A)/(2D/ 

^6=  ^Vp2-'/P|)e3i*  a7=  (3/p22-i/Pi1)e4i, 

A8=I2+5A2/  A9=I2-A2  A|0  = -2lE4/Pl(, 

B, ,  = 12^9-I)/{2(|+A)[(,  + Aj2  + I2j^ 

B2Z-  {(  a0a?  -t  4A3;-  4a(a-u(2o-a;24  I2;]}  I2£4(A-0  2\2!7  (fi-0/pnl 

HPR=(h/Ri)2(50)2//2,  u)2=(2(/<e;2A3-^e;2M4+^^^ 

F,  = K0  ( A4+  D A6  Vc  2 A,) , 

F2  = K0(l  + l>j((A5  tW(A7  + A5)i  A|0  - 42(  A4  + iJA6)/(2^,)JP(  50/e/(  , 

F}  = - (ity)  A2R,/(2Ai$0) 

Tn,  = WI,  >7,  = T7/S0, 

(4,  - 3o(  4 ^,2+  An,  ■+  l )/z,  Ef2  = >*,2/(A2  + m?) } 

3^f/(A2+9mf;y  Er4  = (q-2E2/(4A); 

CT?=  £A/(A2  + 9m?)-A/(A2+mf;3E3;  CTe=  2 ( <r3  - Cr2;  E3 
(^7  = A (l-A)  CEt2  E2  / 2,  <4g  = A mf(  2A2L  nn(z)E2/ f4A(A2+mf  ;], 
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$ 9 "Ari|  ( (q-34  £3^  Gr|0  - C(  I + A ~\2)  (j5  + i>  (qt,  ] yj,  f G[J 

$n  = K9  ( $5  + y $€>)  , $, 2 = 2(Yl2^4i  Cj7  -f  Q-s  ), 

$i3~  ,<f4/  <T/4  = f ™?+'>1?')2B5/(2-(*'iXCi-A)J-+rvi?)3} 

Qr/5=4  5a(m12--*-  An*}^  Grl6  = 4n,m,  ( i + y)  ( E}- 2.)/(  i-p)  , 

Cfl7=2  [ 3h,2/a  + ( £3  -2);  Cr/8  = - n,(»i?+rj*)  E$/(4- 8A), 

<7/9  = O,  (3»r  + d-2D)™iL)')(E1-2)/(i2*-4)i 

F4=  C<t„<2<7,<^- (Tli)-^,o(2Crl^f3-(rn(r/0)]/[2<T,(2(rl(7a’$foJl 

F5-  — < 2 (jr,  {7,3  - Cfn  Crlt> )/ ( 2 - Cj-20  ) ( 

Fg  = < ^/9  <7/6  ~ CT/ g $-,7  )/(<rlSCrl7'  Q-H  ) , 

r9  “ 1 <7/6  <7/8  " <T/9)/ ( <f/f  <7,7  " <4/6  ) , 

F'o  = ("*,(P»?-*nl2-)]  Es  / ( U$e  >nf  (2A-D]j 

£ 

~ ^ 2 n<e>)  ^4  + ^'4"  2<T4^,F5  K8  - ( Grs  + p (r  Kd  FA](SQ)2/e> ,,  t 
V.y-  Ch/Rt)2(  9/n,4-f  9«,vf  2mfnz)Bl0  (56)* / i 156  (+- *)  &,  J; 

^2  = (h/R,)2  ($/$  Fg  -f  Q-/8  F9)  (Sd)4/B,it  j 

^3=  (h/R^htn* -"?)*£*;  fse;4/r/2g^0(^-A)23< j, 

[Qw)it  = ^2 T,  / u5^t 

Qj=  C(0„sej/rR,T,  cost),  f0r  ramp-type 

0,3  — ( ( Q w )$  / ( /^,  i a Co5  T )/  for  sinusoidal  type 

Q't  = F,  Qj  t F3°<  S0/R,  / 

$20~  21  [ ( A - 1/2  )2t  I2J"';  Q-J(=  - 2y  Gt20  q;  £6j 

$21  — 1[(3A-I  )2-H2j;'  ^3=  K0Q3  = 2T(x\l2)> 

$25  ~ ' ^ C 2 u <Tzq  X + A ] E ( 1 ' / $2t>  ~ " £ + 4 ( *v,-f  I J J 

$2-1  = ' C $z^  2a[vlEb  iii+*»)Crlbui-»»crzllt 

$X<=)  ~ + (l<2m,)Z]}  -t  (1-2  (1-2  IV,  l*]'1 

(j3c  = K0  Q3  Gf2^ E?  , ^?,«iUlt<i+ro()2r',  <^2  = lUui-™t)2]'\ 

$33  =0(('‘A)(^  A-d-AjOim))  ^,-KA-d-AUi'iv))  q-32]i'L  2\(a\v2)'1}  E2  , 
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Q'34  = C(A-l/2)2+Cl+m1)2]'1,  Gr35-=  C(A-'/2j2+(  m.-I)2]'1 

= (3A-l){C(3A-i)2+(I  + 2m,')2]"1-  ( C3A-I)2+(I-  2m,)2J~lJ^ 

£-T36  = ~ & i Efe  { C (m  + l)l2-  *>%■]  Cj-^  t ((m-i)I2+^  ^4  ] <q~3$)  , 

6j-3g=  K8Q3  Gt?7ETj 

q-39=  -«C(i-A)/lE2]{Cl(^  + ')  -t  (i-X)A  i"'J^3l  + 

-l-  [l(  w,-l)-A  (i-A)l'1;  <r3Z  - 2w,J  ( Az+  mf)'1 

<r^o  = -9,'Efe  t>+^2JCr2o  ' ^,  = -«  [(i'\)/ie2]{ + 

^42=  2l[C2A'\)24  4(  I + m,)2]*1,  (r43  = 21  [(2,\-i  )2  + 4 ( 

^■44=  Q[&6{((i4m)-i,2;)  ^41  + Cd-^)-y2]  ^43}, 

(445  = CA[fi-A)/l£2]{t"f<+iJ^,4m  ^/z  4 ^Aj  (j]r3  / + 

-t  C-(»+«)n-»»)/2  tyA]  ^32  - 2M/(j\l+m,2)}~\ 

Tc.r  — “ •^■mn  B,,  { ( ( Cjii  - ^Gf25-t  Gr25-)(  A2-+I^if)  4 ( £?2  8 - ^50  + GfjjX  4 

4(Gr3fc  -U  (qr3g  +^)Am,+  C ( <*>0 - Cr2i  + fy)  + (Q44  "fyo  + <T45)/2]»f  JtW*)]' 
Q33  = ~ ( F2/odlU*/+')C  2- 2 CoS  ( C0t,/ 1 )],/2  « 50//?,  , 

Qn  = K9  Qj?  ( A4  + i>A6)/f2iil); 

£t4*>  = _2  w^i«  $/»  ; 641  - Kd  Qh  Q-22  Ej , Q-4i  - a,,  Cr2i/a[ , 

0*4^  ~ Crjb/Q.  , ^50  - K9  q33  £r2?  e/,  Er^i  =■  i<g  Qy5i  Et3j  e7 

^(sz  *3 II  £4-40 /& I , £^53  “ £5  11  £t44  /Q\ , 

2 

~£V  B,,  f ( ( - v Cf^j  )(\  -t  mf)  4 ~l>  Cffo)  (m?-J\x)/2  4 

( *444  - o Cr5|)A  m,  + i(  Cf5z~  <T5i~  (jro)/zJ  nf  J (SO)2}  . 


APPENDIX  II 


FORMULATIONS  OF  THE  CYLINDRICAL  SHELL 


~ vL/rrrjiR,  F7  = ~4-(  t + D)R/T\, 

COc  = i-^2+Ch/R)2(R/L)47rV/2  + (SV2/^2), 

7r?  = W7 T/l, 

HRP  = (h/R)2  (ynzR2  + n2)2/tZ  , 

A = u (?Pj2l2-S)R/('R)L2)i 
A,  = >7  2 + a^2R2 

J\-Z~  i('+'J)/2v')'nAI 

A3-  ™2Rz- (8R2/lz)  i 5U-u)-h2/e  - [ ( 8(i'L»»2)/(fr2L2)], 
A4  = 4 *n2F(2i  4An2 


A 5=2  n + u;  t-iamR 
Afe  =[?r>R  A + 3nzJ4, 


A 7 = tti R f ?n2Rz-t  ri2)/4y 
Ag=  ??  C 3t?2 ^2R2J/4, 


RR  4A,-)iA2,  B6  — A3A,- CC  = >7  A 3 - A A£; 
(n^n)2=  i I •+  HRP  - ( >7  CC  +ARfl)/B6], 

= (ll/R)2  (8  >ri4R4  t 4 n4-  1 2 >v2 r>2 Rz ) / 32. 

ylic.~  2-(h/R)  ( A7  A(,~  2 AfA-j  A q -t  A4  Aq)/(  a4  At,  - A$), 
= ( h/R)2  ( P"2-  wVj2//6^ 

<T6o  = - [«(t+j,)-(2U/Tr)F7c<Afl(Co‘-.nf))]  ^2R2 


for  sinusoidal 
r load 


^61  = ~ C«  C * + L» ) - 2 [(  4m2-v2)/(4rr)2-i)]F7o</C7jR  (A?-Sl2)]')  ?i2 

Q"fc°  — “ ^ (z  u/t[)  F-^  / ( R cA2  )]  >7i  2 R2  1 for  ramp-type 

Q'foi  = (i-tD)  - 2[(4»n1-^2/)/<'4ni -i;J  F7o(/(jrRcC)(]]  772 

* / 
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Tccr  =-[(il^)2/(Grfoo-fGrfo,)j 

= ( 2 p/'tt)  q o(  (i/cdct0')(2-2  CoS  ( cJcto/t )]  1/2 
Cf64=  »ZV ( Cr63/u)(n2/(fn'2F\2)); 

— ( ( -^-ron)  /(  3 + Q"64  )3  • 


2R1/(R  ode), 


APPENDIX  III 


FLOW  DIAGRAM  OF  COMPUTER  COMPUTATIONS 


CONICAL  SHELLS 


SINUSOIDAL  TEMPERATURE 


RAMP-TYPE  TEMPERATURE 


FOR  CYLINDRICAL  SHELLS,  THE  PROCEDURE  IS  THE  SAME  EXCEPT 
R/h  AND  ARE  REPLACED  BY  R/l  , R/h  AND  u)c  AND  THE  DO 
STATEMENT  FOR  0 IS  DROPPED. 
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